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Abstract. By means of the Density Matrix Renormalization Group technique, we have studied the region 
where XXZ-like behavior is most likely to emerge within the phase diagram of the F-AF anisotropic 
extended (J — J') Heisenberg chain. We have analyzed, in great detail, the equal-time two-spin correla- 
tion functions, both in- and out-of- plane, as functions of the distance (and momentum). Then, we have 
extracted, through an accurate fitting procedure, the exponents of the asymptotic power-law decay of the 
spatial correlations. We have used the exact solution of XX Z model (J' = 0) to benchmark our results, 
which clearly show the expected agreement. A critical value of J' has been found where the relevant 
power-law decay exponent is independent of the in-plane nearest-neighbor coupling. 
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1 Introduction 

The spin- 1/2 extended Heisenberg model is one of the 
most studied spin systems in the literature as it is the 
prototype of a geometrically frustrated quantum system 
and allows to study the interplay between two funda- 
mental aspects of spin physics: quantization/fluctuation 
and frustration. Such competition is often responsible for 
the appearance of exotic phases and anomalous phenom- 
ena whose explanation is still under debate. The one- 
dimensional realization of the extended Heisenberg model 
(hereafter, J stands for the nearest-neighbor exchange cou- 
pling and J' for the next-nearest-neighbor one) has been 
mainly studied in the AF-AF regime (J > 0, J' > 0) 
according to the many experimental realizations. Only re- 
cently, the F-AF regime ( J < 0, J' > 0) has come into the 
limelight as it has been proposed to explain the anomalous 
experimental behavior of compounds containing chains 
formed by edge-sharing CuOi plaquettes: Rb2Cu2Mo30i2 
[U, NaCu^O^ |2], and LiCuVO^^ In such class of com- 
pounds (both edge- and corner- sharing ones), the sign 
and the absolute value of J depend sensitively on the Cu- 
O-Cu bond angle and on the distance between copper and 
oxygen ions. In the compounds mentioned above, the Cu- 
O-Cu bond angle is slightly larger than 90°. Accordingly, 
the usual antiferromagnetic nearest-neighbor exchange be- 
tween Cu ions becomes ferromagnetic (the 02p orbital 
hybridizing with the 3d orbital of one Cu ion is almost 
orthogonal to that of the next Cu ion). This nearly or- 
thogonality makes the nearest-neighbor coupling in the 
edge-sharing case smaller, by more than an order of mag- 
nitude, than in the corner-sharing case. In edge-sharing 
compounds, the next-nearest-neighbor exchange couphng, 
through the Cu-O-O-Cu path, is generally antiferromag- 



netic, usually with a magnitude of a few tens Kelvin. 
Despite its smallness, the next-nearest-neighbor exchange 
coupling interaction has a pronounced effect on the physi- 
cal properties of these systems since the nearest-neighbor 
exchange coupling is also very small owing to the nearly 
orthogonality {J' / J ~ —1/3). On the theoretical level, 
it is known that for J'/ J > —1/4, the ground state of 
the system is ferromagnetic 4^. At exactly J'/ J = —1/4, 
a new singlet ground state, degenerate with the ferro- 
magnetic one, appears and has a resonating-valence-bond 
form (under periodic boundary conditions: Ref. [5]). For 
J' / J < —1/4, only speculations are available, which range 
from a gapless singlet ground state to a gapped one 



In this manuscript, we have studied the anisotropic 
version of the spin- 1/2 extended Heisenberg chain with the 
intent of locating the region in its phase diagram where a 
XXZ-like behavior can be observed despite the presence 
of a finite next-nearest-neighbor coupling J'. In particular, 
we have analyzed the dependence of the critical exponents 
of the in-plane and out-of-plane correlation functions on 
J'. Quite unexpectedly, a critical value of J' has been 
identified where the in-plane critical exponent does not 
depend on the in-plane nearest-neighbor coupling and the 
out-of-plane correlations qualitatively change their behav- 



2 



Adolfo Avella et al.: XXZ-like phase in the F-AF anisotropic Heisenberg chain 



2 Definitions and Method 



The anisotropic version of the F-AF Heisenberg chain reads 
as 

i i 

+ j'Es,s,+2 (1) 

i 

where the nearest-neighbor out-of-plane couphng constant 
Jz has been taken ferromagnetic (Jz > 0) and the next- 
nearest-neighbor couphng constant J' has been taken anti- 
ferromagnetic (J' > 0). Results can be easily mapped be- 
tween negative and positive signs of the nearest-neighbor 
in-plane coupling constant J± that has been fixed positive 
(antiferromagnetic). Hereafter, will be used as energy 
unit {Jz — 1). The Hamiltonian model ([T]) opens up the 
possibility to analyze, in great detail, the role/relevance 
of the geometrical frustration in enhancing/quenching the 
quantum fluctuations. In fact, it contains both J' =/= 0, 
which can drive either gapless liquid [51|7] or gapped dimer- 
ized phases and Jz ^ J±, which can obviously drive 
quite different responses in- and out-of- plane. The very 
rich phase diagram, the competition between frustration 
and anisotropy and the absence of evident small param- 
eters has attracted the attention of many researchers in 
the field. When Jz, J± and J' are all of the same order 
of magnitude, exact (or quasi-exact) treatments, such as 
Bosonization and Renormalization Group, are inapplica- 
ble and numerical techniques are the only possible tool 
of investigation [9]. In this manuscript, we have used the 
Density Matrix Renormalization Group (DMRG) tech- 
nique ^lOllllj in order to obtain the ground state properties 
of the Hamiltonian ([T|). We have studied a, L — 100-site 
chain, subject to open-boundary conditions, and we have 
retained up to the lowest 200 eigenstates of the reduced 
density matrix in the basis of each DMRG block. The 
real-space spin-spin correlation functions [both in-plane 
{{S^'^Sy}) and out-of-plane {{S^S^,})] have been calcu- 
lated, as functions of the distance d = \r — r'\, by choosing 
r and r' symmetrically with respect to the center of the 
system. The momentum-space correlation functions have 
been obtained as the Fourier transform of the real-space 
ones. It is worth noting that Fourier transform is gen- 
erally useless when open-boundary conditions are used: 
functions in direct space generally depend on two set of 
coordinates and the transformed functions generally de- 
pends on two momenta. However, larger is the system 
under analysis, smaller are the non-diagonal terms in the 
transformed function (they are exactly zero when periodic 
boundary conditions are used) . We have carefully checked 
that our data satisfy this condition: non-diagonal terms in 
momentum space are safely negligible in the whole region 
of parameter space under analysis. Such test has confirmed 
that 100 sites are sufficient to our purposes and that the 
effects of open boundary conditions are under control. 
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Fig. 1. XXZ-like region in the (J', J_i) plane. The solid line 
is a guide to the eye. The dashed line is the classical transi- 
tion line between an in-plane antiferromagnet and an in-plane 
incommensurate spiral phase. 

3 Results 

In the region under analysis {J±/Jz > 1, J'/Jz < 1), the 
classical phase diagram (i.e. the one obtained for classical 
vectors of constant length) presents a transition between 
an in-plane antiferromagnet (main momentum is fc = tt) 
and an in-plane incommensurate spiral phase with main 
momenta k — tt zt arccos( Jj_/4 J') for J' > Jl/4 (dashed 
line in Fig.[T]). Quantum fiuctuations modify this picture. 
Quasi-long-range correlations, with power-law decay at 
large distances, replace the long-range correlations with 
constant values. In the simple J' = case {XXZ model), 
as soon as the rotational invariance is broken, the long- 
distance correlation functions assume a simple power-law 
form [12] : 

(S'oS'r) {-lyAzr'/^' + BzT-^ 

where rj — — &vccos{J z / J i) / , 1 < J^jJz, has been 
determined in Ref. |13j . 

We have analyzed the equal-time two-spin correlation 
functions, both in-plane {S^S^) and out-of-plane {SqS^.), 
as functions of the distance r, together with their Fourier 
transforms, in order to easily recognize short-, quasi-long- 
and long- range orders and their characteristic features. A 
long-range order is characterized by a correlation in mo- 
mentum space showing a Dirac delta (a peak of height 
proportional to the size of the system) at some momen- 
tum k, which corresponds to a correlation in direct space 
with constant amplitude and modulation e''^''. The system 
under analysis shows such a behavior only in the ferro- 
magnetic phase {J±/Jz < 1, J'/Jz ^ 0.25). A quasi-long- 
range order is characterized by a correlation in momentum 
space showing a power-law divergence 1/ |A: — fc|'' (a peak 
on a finite system) with exponent ^ 77 < 1 at some mo- 
mentum k. Such a behavior corresponds to a correlation 
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Fig. 2. Equal-time two-spin correlation functions in-plane {S^S^) (left panels) and out-of-plane {S^S^) (right panels) as 
functions of the distance r (top panels) and momentum k (bottom panels) for J± = 1.2 and different values of J'. 



in direct space showing a power-law decay 1/ |r| with 
exponent < 1 — 77 ^ 1 and modulation e''^''. We have 
extracted k from the correlations functions in momentum 
space and rj from those in direct space. In fact, the re- 
lation l/|fc|'' ^ holds only for infinite systems 
and is badly violated for finite systems. A short-range or- 
der is characterized by a correlation in momentum space 
showing a Laurentian at some momentum k of width A, 
which corresponds to a correlation in direct space showing 
an exponential decay e~^^ with characteristic length 1/A 
and modulation e''^'". 

In Fig. [2] we have reported the equal-time two-spin 
correlation functions, both in-plane (left panels) 

and out-of-plane (S^S^) (right panels), as functions of 
the distance r (top panels) and momentum k (bottom 
panels) for J± — 1.2 and different values of J'; results 
differ only quantitatively for higher values of J±. Up to 
J' = 0.32, the in-plane correlations clearly show a simple 
quasi-long-range (linear in log-log scale) antiferromagnetic 
order (left panels): a well defined peak centered at fc = tt 
is clearly visible. At J' = 0.34, the in-plane correlations 
lose the quasi-long-range character: the decay is exponen- 



tial (top left panel and its inset: linear in semi-log scale) 
and the peak at A: = tt becomes a Laurentian (bottom left 
panel). The out-of-plane correlations (right panels), com- 
pletely absent in the classical case, clearly show a linear 
well (|fc| — > r~^) at J' = that changes its concavity on 
increasing J' (bottom right panel) . A small peak develops 
at fc = ±7r for 0.2 ^ J' < 0.34. For larger values of J', 
the out-of-plane correlations develop a small (short-range 
— linear in semi-log scale) antiferromagnetic component. 
If we use the appearance of an exponential decay in the 
in-plane correlation functions as the signature of a tran- 
sition to a different phase, we can locate the XXZ-\ike 
region in the (J', J±) plane and draw the transition line 
shown in Fig. IT]) (circles and solid line). It is worth not- 
ing that, for the largest reported values of J' (J' ^ 0.30), 
the results for the out-of-plane correlations in direct space 
fall on three different curves (top right panel) and, corre- 
spondingly, in momentum space (bottom right panel), it 
is clearly visible the emergence of peaks in proximity of 
k = ±7r/2. Such behavior is most probably connected with 
a spiral phase, reminiscent someway of the classical one, 
appearing beyond the the massive phase that we are using 
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Fig. 3. Exponent r] of the power-law decay over distance of (S''~S''~) for different values of J± (left panel) and J' (right panel). 
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to bound the XXZ-like phase under analysis. We will not 
discuss further neither the massive nor the spiral phase as 
they are out of the focus of this manuscript. 

It is worth noticing that the observed behaviors are 
quite well described by (pi. Then, if we do believe that the 
Q do not change for small enough values of J' , our analy- 
sis simply accounts for the dependence of the critical expo- 
nent T] on the value of J'. While the in-plane correlations 
have as leading term the antiferromagnetic one for any fi- 
nite value of 77, the out-of-plane correlations present a com- 
petition between the two components (antiferromagnetic 
and ferromagnetic) depending on the value of rj. In par- 
ticular, for values of 1] larger than —1/2 (—1/2 < 77 ^ 0), 
the ferromagnetic component is dominant, while the anti- 
ferromagnetic one prevails for 77 < —1/2. Now, for J' — 0, 
77 can only decrease from to — 1/2 on increasing J± . Ac- 
cordingly, the out-of-plane correlations can acquire an an- 
tiferromagnetic component that has an overall exponent 
{l/rj ^ —2), at most, as small as the ferromagnetic one 
(—2). For J' > 0, 77 can reach (see Fig.[3|, for a fixed value 
of J_L, much lower values than —1/2 and, in principle, it 
becomes possible to explore cases in which, for the out- 
of-plane correlations, the antiferromagnetic component is 
dominant. 

In Fig. |3] we have reported the critical exponent 77 of 
the power-law decay of the in-plane correlation functions 
for different values of J± (left panel) and J' (right panel) . 
In the relevant parameter region, the exponent 77 has been 
extracted through an accurate fitting procedure of the in- 
plane correlation functions in direct space (in log-log scale 
plots). In order to benchmark such a fitting procedure, we 
have also analyzed the exactly solvable J' = case and 
compared our results with the exact ones, resulting in an 
excellent agreement (right panel). It is clearly visible that 
there exists a critical value J'^ ^ 0.3 where 77 is indepen- 
dent of J±. It is worth noting that for J' — J'^, -q is very 
close to —1/2, the critical value for the out-of-plane cor- 
relations (ferromagnetic vs. antiferromagnetic). 

In Fig. [§ we have reported the behavior of both in- 
plane (left panels) and out-of-plane (right panels) correla- 
tion functions in momentum space for values of J' below 
(J' = 0.2) and above (J' = 0.3) the critical one in order 
to extract any significant difference in the behavior across 
such a transition. Below J^, the in-plane correlation func- 
tions show a monotonic decrease of the peak intensity on 
increasing J±. Below J^, the peak generally has a lower 
intensity with respect to the former case. Moreover, and 
more interestingly, the intensity of the peak first increases, 
on increasing J± , and then decreases and rapidly saturates 
to values similar to those obtained below J^. The out-of- 
plane correlation functions simply show, at first glance, 
a systematic increase of their value at ±7r on increasing 
J±, for both values of J'. What is really reminiscent of 
the behavior found for (S-^S-^) is instead the slope of 
(S^S^) at ±7r. This latter changes sign at exactly the same 
value where the maximum of the intensity is reached for 
(S-^S-^). Moreover, for those values of J± where the in- 
tensity of the peak in (S-^S-^) increases, two well-defined 
shoulders appear in (S^S^). 



4 Conclusions 

In conclusions, we have found a region in the phase di- 
agram of the spin-1/2 extended anisotropic Heisenberg 
chain where a XXZ-like behavior persists. We have ex- 
tracted the exponent of the power-law decay and found a 
critical value of J' where the former is practically inde- 
pendent of J±. 
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